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I. INTRODUCTION 

Since the Hawking effect^ on black hole was found, the evaporation of black hole 
has been investigated in several coordinates by miscellaneous methods such as path 
integral approach,^ tortoise coordinate (r*) transformation,^'^ and r*-coordinatc ana- 
lytical extension,^ etc. Among these methods, the generalized tortoise transformation 
method has been used widely in the discussions not only on evaporation of static black 
hole and stationary black hole, but also on that of non-stationary ones.^ Much more 
progress has been made. But this method can not give an exact solution of the radial 
(r*) equation, the radial wave function, which can be analyzed only in asymptotic 
expression. 

Couch^ obtained a series of exact solutions by transforming the separated radial 
equation into modified Whittaker-Hill equation under some special conditions. But 
these solutions seem to have nothing to do with the discussion on black hole evapora- 
tion. 

Solutions to generalized spheroidal wave equation have been studied to some extent^' 
by using power series expansions around regular singular points, so that three-term re- 
currence relations between coefficients can be manipulated in terms of the continued 
fraction method. Leaver^ has shown that Teukolsky master equations in Kerr geometry 
are, in fact, spin- weighted generalized spheroidal wave equations. 

It appears to be more important to obtain an exact solution to the radial equation 
for this is crucial in discussing Hawking effect of black hole. However, it is very difficult 
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to do so. It is this motivation that stimulates our present research. The main aim of 
this paper is to show that the separated radial part of a massive covariant Klein-Gordon 
equation on the Kerr-Newman black hole ( KNBH ) background is a generalized spin- 
weighted spheroidal wave equation of imaginary number order. In this paper, we shall 
discuss the solutions to a massive complex scalar field in the KNBH geometry with 
three parameters. In the non-extreme case, its general solutions of the separated parts 
are spin-weighted generalized spheroidal wave functions®^® and some special solutions 
to the radial equation with physical interest are given. General solutions to the radial 
equation in the extreme case shall be briefly discussed. Finally, we show that the radial 
equation in the extreme case is a confluent equation of that in the non-extreme case. 

Sec.II deals with variable separation of a sourceless complex scalar field on KNBH 
and solutions to the angular part. In Sec. Ill and IV, the radial equation is solved 
in both non-extreme and extreme cases respectively. In Ill(a) we reduce the radial 
equation to standard form, and in Ill(b) and (c) we obtain general solutions and 
special ones including case {cu — /i — 0) respectively. Five-term recurrences between 
coefficients of solutions in power series forms are given in both cases. In addition, 
we give solutions in integral forms and some special solutions of physical interest in 
the non-extreme case. Conditions for general solutions exist are given in these two 
cases. Sec.V is devoted to discussing the connection between the radial equation in the 
extreme case and that in the non-extreme case. Finally, we point out some probable 
applications and generalization of exact solutions in Sec. VI. 
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In appendix part, three-term recurrence relation between coefficients in power series 
expansions around regular singular points for generalized spheroidal wave equation are 
presented. 



II. SEPARATION OF KLEIN-GORDON EQUATION AND SOLU- 
TION TO THE ANGULAR EQUATION 



The Kerr-Newman line element and electromagnetic one-form are given in the 
Boycr-Lindquist coordinates as foUows^'^ 

ds^ — g^j,jydx^dx^ 



with event horizon function A = — 2Mr + + e^, and E = + cos^ 9, where 
mass M, charge e, specific angular momentum a — J/M being three parameters to 
describe KNBH. ( Use Planck units system G — h — c — 1, and denote — ). 

The determinant of KNBH metric tensor is g = dei^g^y) = — E^sin^^, while the 
electromagnetic four-vector potential apparently satisfies the following covariant 
Lorentz gauge condition: 



In curved spacetime, a sourceless scalar field $ with mass // and charge q obeys the 
covariant Klein-Gordon equation ( KGE ): 



A 






1 



(3) 



(□,-/i^)$ = 0, 



(4) 
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where d' Alembert operator on KNBH background is given by 

= ^{^[(r' + a')dt + ad^ + iqerf + a,(Aa,) 

+(asin^9t + ^9^)' + -^9,(sine9,)}, (5) 

here covariant gauge differential operator being — d^ — iqA^. 

The scalar wave function $ for KGE of Eq. (4) has a solution of variables separable 
form r, 9, (p) = R{r)S{e)e'^""^-'^*'^:^ 

■^[uj{r^ + a^) - qer - ma]^$ + dr{Adr^) - /i^E$ 

771 1 

-(aa;sin^ + -deismOde^) = 0. (6) 

sin 9 sin 6/ 

The separated results of the above equation are 



. [a;(r^ + a^) — qer — ma] 
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dr[AdrR(r)] + { '^'' ' - ^ _ ^2(^2 ^ ^2^ _ ^ ^ 2mau;}R{r) = (7) 



and 



. ^ae[sin^ae5(^)] + [A-^^ + (/x'-a;2)^2g.^2^|^^^^^Q ^g^, 
smc^ sm 9 

where A is a separation constant. 

The general solutions to the angular part are ordinary spheroidal angular wave 
functions^'^° with spin- weight s — 0. When a'^{uJ^ — jj?') — 0, these solutions degenerate 
to Legendre spherical functions. 

Let X = cos^,5(^) = S{x) = (1 - x2)'"/2e(a;), Eq.(8) should take the following 
forms: 

2 

(1 - x'')S"ix) - 2xS'ix) + [A - + aHu"" - /i^Mx^ - l)]S(x) = (9) 

1 — x^ 



and 

(1 - x^)e"{x) - 2(1 + m)xQ'{x) + [A - m{m + 1) + a^{u;^ - f/) (x^ - l)]e(x) = 0. (10) 

Here and after, S'{x) — dS{x)/dx, etc.. 

The eigenfunctions to Eqs.(9) and (10) are generalized spheroidal wave functions^"^ 
S{x) = S'7'°(c, x) with eigenvalue A = + = a^(//^ -a;^). When = 0, Eqs.(9) 
and (10) are special cases (s = 0) of the following spin-weighted spheroidal wave 
equations^~^'^° 

(1 - x'')P"{x) - 2xP'{x) + [a^u^x^ - 2aujsx - ^ ^"^^^ + s + \']P{x) = (11) 
and 

(l-x^)g"(a;)-2[s+(l+m)a;]Q'(a;) + [A'-(m-s)(m+s + l)-2acusx+a^tuV]g(x) = 

(12) 

where P{x) = (1 - x)l'"+^l/2(l + x)\'^-'\I^Q{x) and x = cos^. 

When ao; = 0, the solutions to the above equations are Jacobi ultra-sphere D- 
functions^°£)^^(a;) or spin- weighted spherical harmonic functions^^ with eigenvalue 
A' = l{l + 1) — s{s I — max(|m|, |s|). In general case, the solutions should be the 
generalized spin-weighted spheroidal wave functions®"® P(a;) = P^^g{c,x), = —c^LtP'. 
By taking account of some reasonable boundary conditions, these solutions could be a 
set of orthogonal polynomials. 

In the following, we shall assume that all parameters, M, e, a, /x, g, m, are nonzero, 
and discuss the radial equation of Eq.(7) according to two cases, namely the non- 
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extreme case (M^ 7^ + e^) and the extreme case (M^ = + e^). Special case 
(a; = // = 0) will be included in subsection III(c). 

III. SOLUTIONS TO THE RADIAL EQUATION IN THE NON-EXTREME 

CASE (M2 ^ + e^) 

(a) Simplification of the radial equation in the case {e 0) 
In this case, we put e — \/M'^ — — e^, (0 < £ < M). After making substitutions 
oir = M + ez and R{r) = R{z) = {z - 1)^1^+^1/2 ^ iy\^-^\/^F{z), the exterior 
horizon and interior horizon are located at points r± — M ± e, {z — ±1), respectively, 
the radial equation of Eq.(7) can be reduced to the following standard forms: 

{z^ - 1)R" + 2zR' + [e^co^ - n^){z^ - 1) + 2e{Au; - MpL^)z + 

+{2u;^ - H^){2M^ - e^) - 2qeMuj - X]R ^ {1< z < 00) (13) 

and 

(z^ - 1)F" + 2[iA + (1 + iB)z]F' + [e^{uj'^ - i/){z'^ - 1) + 2e{Auj - Mi/)z 

+ (2a;2 - i/){2M^ - e^) - 2qeMu; -X + A^-B^ + iB]F = (1 < z < 00) (14) 

where A = 2Muj - qe, sB = u;{2M'^ - e^) - qeM - ma. 

In order to study behaviors of solutions to Eqs.(13) and (14) in the interval (—1 < 
z < 1), we rotate firstly T from real axis to imaginary axis T — it after making 
substitution oi z = coshT = cosh(iT) = cost, then return to real z-axis z = cost. 
Therefore, Eqs.(13) and (14) have corresponding forms in the interval (|2;| < 1) as 
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follows 



(1 - z^)R' 



ill 



2zE! + [e\u^ - i/){z^ - 1) + 2e{Auj - Mii^)z 



{Az + Bf 

1-^2 



l/){2M^ - e^) - 2qeMcu - X]R = 



{-l<z< 1) (15) 



and 



(1 - z'')G" - 2[A + (1 + B)z]G' + [s\uj^ - i/){z' - 1) + 2e{Au - Mf/)z 

+{2ou^ - n^){2M^ - e^) - 2qeMu - X + - - B]G ^ {\z\ < 1) (16) 



The essence of our manipulation is that we extend the domain of z from the real 
axis to complex 2;-plane z — x + iy ai first and then make an analytical extension on 
complex 2;-plane from region outside unit circle {\z\ > 1) to that inside it {\z\ < 1). 
The handled issue is that only derivative terms change a negative sign, while the non- 
derivative terms, namely terms in square brackets, make no change in symbol. This 
method is equivalent to that Eqs. (13-16) are solved initially and then the solutions are 
made analytical extension on the complex z-plane. 

Both Eqs. (13) and (14) are generalized spin- weighted spheroidal wave equation^"® 
with imaginary number order, while both Eqs. (15) and (16) with real number order. 
The formers are suitable especially to study problems about scattering state, whereas 
the latters more convenient to investigate energy levels of bound states. Furthermore, 



where we have made a function transformation 



R{z) = il-z) 



B + A 
2 



{i+zy-^G{z). 
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we can apparently find out connection between poles of scattering amplitudes and 
energy levels of bound states. Actually, domain in which z takes values in Eqs.(13) 
and (14) is on the x-axis, while those in Eqs.(15) and (16) on y-axis. So these equations 
can be thought of equivalence. However, to be convenient, we have made a restriction 
on intervals that z takes values of 1 < 2; < oo in Eqs.(13) and (14), while of \z\ < 1 in 
Eqs.(15) and (16). 

When /i = or M = 0, if taking R\{z\) as the first solution to Eq.(15) in the 
interval of \z\ < 1, then R2{z2), the second one to the same equation in that of > 1, 
might be: 

+1 

R2{z2) = {Z2 - l)^(-22 + 1)^ / e™^^(l - -2i)^(l + z,)^R{z,)dz^. (17) 

-1 

Here, we have assumed that a; > 0. Integral equation of Eq.(17) connects irregular 
solution R2{z2) with regular solution Ri{zi). 

Comparing Eqs.(9-12) with Eqs.(13-16), especially Eqs.(ll, 12) with Eqs.(15, 16), 
we can draw a conclusion that the separated angular and radial equations are ordi- 
nary differential equations of the same type, generalized spin-weighted spheroidal wave 
equations. Furthermore, we discover that e, A, B correspond to —a, — s, —m respec- 
tively when = 0. There may exist three pairs of power series solutions to generalized 
spheroidal wave equations around singular points z — ±1, 00 respectively. Added with 
some proper boundary conditions, these power series expansions of spheroidal wave 
functions can be cut off to be polynomials. 

Therefore, in the following subsection, we shall only study the generalized spheroidal 
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wave equation. The reader who has more interest in this equation, can find more in- 
formation in Refs. 6 and 7 ( and references cited therein ). 
(b) General solutions to the radial equation {e ^ 0) 
The standard generahzed spin-weighted spheroidal wave equation that we reduce 
to study is as follows 



(1 - z^)W"{z) - 2[a + + l)z]W'{z) + [-f^{z'^ -l) + 2Sz + X- P]W{z) = 0, (18) 



where A is a redefined eigenvalue which could make W{z) finite at z = ±1, and region 
of taking values could be the whole complex z-plane. 
(i) For the radial equation of Eq.(16), we have 



(ni) For the angular equation of Eq.(lO), we have 

a = 0, /3 = m, 7^ = a^(u;^ - /x^), 5 = 0. 

The form of Eq.(18) is invariant both under Laplace-transformation and by chang- 
ing parameters a, /3, 7^, 5, z into —a, /3, 7^, —5, —z respectively. Namely, W{z) — W{a, /3, 7, 
satisfies the following integral equation 



a = A, /? = 7^ = £2^a;2 _^2^^^^ ^^^^ _ ^^2-) . 



(n) For the angular equation of Eq.(12), we have 



/3 2 2 2 c- 

a — s,p — m,^ —auj,d — — 



ascu; 








W{ 



-6 



■/3,7,q;7; — ), 



(19) 



7 



10 



W{a,(3,j,d;z) = W{-a, (3,j, -S; -z) 
+00 - 
^ f e-'^'W{ — ,-(3,j,ar,t)dt. (20) 

i ^ 

The above formulae are integral solutions to Eq.(18). If one knows a solution, then 
he can obtain another by integral transformations of Eqs.(19) and (20). It is obvious 
that solutions are symmetry under the following condition: 

a = 5 = 0, (7 7^ 0). 

This just is the case (iii). At this moment, the symmetric solutions are ordinary 
spheroidal angular wave functions. 

Now, we consider a solution to generalized spheroidal equation of Eq.(18) which is 
in power series form in the interval of — 1 < 2; < 1. According to the knowledge of 
ordinary differential equation, one can know that Eq.(18) has two regular singularities 
{z — ±1) and one confluently irregular singular point {z — 00). As z = is its 
ordinary point, we can make a Taylor expansion of W{z) in the vicinity of ordinary 
point {z — 0): 

00 

W{Z) = Wniz) = J2 CinZ'\ {\z\ < 1). (21) 
n=0 

Substituting power series of Eq.(21) into Eq.(18), we obtain five-term recurrence 
relations between coefficients as follows 

(A - 7^ - (3)ao - 2aai + 2a2 = 0, 

25ao + [A - 7^ - /5 - (2 + 2(3)]ai - Aaa2 + 03 = 0, 

7^ao + 25ai + [A - 7^ - /3 - 2(3 + 2/3)]a2 - Gaaa + 12a4 = 0, 
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7^a„_2 + 25a„_i + [A - 7^ - /? - n(n + 1 + 2/3)]a„ 
-2(n + l)Q;an+i + (n + 2)(n + l)a„+2 = 0. 
Redefine coefficients: 

7^ 



A - 72 - /? - n(n + 1 + 2/3) ' 

25 

A - 72 -/3-n(n + 1 + 2/3)' 

-2(n + 1)« 
A - 72 -/?-n(n+ 1 + 2/3)' 
(n + 2)(n + 1) 



A - 72 - /3 - n(n + 1 + 2/3) ■ 
When taking limits n — > 00, we fiave An,Bn,Cn — > 0, and — > — 

Tlien, five-term recurrence relations become 

AnO'n-2 + Bnttn-l + a„ + CnCln+l + Dnan+2 — 0. 



(22) 



After arranging coefficients An, Bn,Cn, and making up them into a quasi-diagonal 
band matrix A and oq, ai, ■ ■ ■ , a„, ■ ■ ■ into a column vector a = (ao, ai, ■ ■ ■ , a„ 
,•••), the above recurrence relations become an infinite tridiagonal matrix equation: 

Aa = ha. (23) 

The condition for solutions of Eq.(23) exist is that determinant of A is zero, 



det(A) = 



1 


Co 




















Bi 


1 


Ci 


Di 














A2 


B2 


1 


C2 


D2 














As 


B3 


1 




Ds 














An 




1 




Dn 






0. 



(24) 
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In fact, this condition could be satisfied, and we have det(A) — > when n ^ oo. 

Matrix equation of Eq.(23), together with determinant equation of Eq.(24) de- 
termines coefficients Oq, Oi, • • • , ■ " eigenvalue A, hence, eigenvalue A will be 
a complicated function of q;,/9, 7,5, as well as n. The second power series solution 
around the same point z — can be obtained by Probenius's method. To be finite at 
z — ±1, power series W{z) could be truncated to be polynomial, and a, 7, S could be 
integers or half-integers. While in general case, solutions to spin-weighted generalized 
spheroidal equation of Eq.(18) are transcendental functions. ^'^ 

Absolutely, solution Wn{z) — Wn{a, (3,j,d; z) of Eq.(18) can be orthonormalized 
to constitute a set of complete functions. 



Solutions Wn{z) at infinity can have asymptotic forms Wn{z) — >• e^^^, {z — >• ±00,7 > 
0). This is consistent with that the Minkowski spacetime is an asymptotic spacetime 
of the Kerr-Newmann black hole. Thus, in-going wave and out-going wave at infinity 
can take form of plane waves. 

(c) Special solutions to the radial equation in the case {e 0) 
In this subsection, we will base our discussion upon Eq.(15), namely 



1 




(25) 



-1 



(1 - z^)R"{z) - 2zR'{z) + [72(^2 - 1) + 2Sz 



{13 + azf 
z^ - 1 



+ (26^2 _ ^2^)(^2M2 - e^) - IqeMuj - }^R{z) = 0, (|^| < 1) 



(26) 



where 
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. ^ , , ^ ^ uj{2M^ - e^) - qeM - ma 

e 

Case-1: when j — S — 0, there exist three situations: 

i) a; = ±/x = qe/M ^0,{a^ 0); 

u) u; — fjL — qe/M — 0,{a — 0) ( This case can be thought as special one in the 
above- head case. ); 

iii) a; = /X = 0, qeM 0, (a 7^ 0). 

Solutions in situations i) and iii) are Jacobi ultra-sphere functions R{z) — p(fi+°''(^-<^) 
(^)/° whereas solutions in situation ii) degenerate to be Legendre functions, R{z) — 
P^{z), or Qi{z). 

Case-2: when MjjF' = 0, a; 7^ 0, 5/{euj) — a, this case has been considered in detail 
by E. W. Leaver.^ 

Case-3: when a = 5 = 0, 7 7^ 0, Eq.(26) is an ordinary spheroidal wave equation, ^'^ 
and its solutions are Prolate spheroidal angular wave functions R{z) = ,Sf'°(7,^). 

Obviously, all these solutions are special cases of general solutions (^^^ S] z) 

= (1 - z)(/3+")/2(i + z)('3-")/2pi/^(a, /5, 7, 5- z). 

Solutions in case-1 will be particular important in physics to black hole evaporation, 
as scattering cross section, stationary state energy levels, emission coefficients of black 
hole radiation, etc., could be analytically computed at exact theoretical level by use of 
Jacobi polynomials. Furthermore, there maybe exist special symmetry in such case. 

IV. SOLUTIONS TO THE RADIAL EQUATION IN THE EXTREME 
CASE (M2 = + e2) 
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In the extreme KNBH case (e = 0), we make substitution r = M{l + x), then event 
horizon is located at a single point (r^ = M), namely x — 0, hence the radial equation 
of Eq.(7) can be transformed into the following confluent equation: 

x'^R"(x) + 2xR'{x) + - ii'^)M^x'^ + 2{Auj - Mii'^)Mx 
+(A + f )2 + (2u;2 _ ^2-) ^2M2 - e^) - 2qeMu - \]R{x) = (27) 

where A = 2Muj - qe, MB ^ Be ^ uj{2M'^ - e^) - qeM - ma. 
Defining 

= M2(a;2 _i/-j^D = M{Auj - Mi/), 
Ae = (2^2 _ ^2^)^2M2 - e^) - 2qeMuj -\-^ + A^. 
and making substitutions: 

then, Eq.(27) is transformed into the generalized Whittaker-Hill equation ( GWHE ) 

- v-^H'\i) + [C2e2i,.€ ^ 2L)e*'^« + 2ABe-*^^ + B^e-^*^^ + Ae]i/(0 = 0. (28) 
Solutions of GWHE of Eq. (28) can be regarded formally as 

+00 

^(0= E ^ne-'^^ n = 0,±l,±2,-... (29) 

n=—oo 

Substituting Eq.(29) into Eq.(28), we obtain five-term recurrence relations between 
coefficients 

^9^-2 + aiP^n-i + (Ae + n'')gn + 2A%„+i + B''gn+2 = 0, 
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+00 n+2 

E'g^hJ^ E ^--"^n, (30) 

n=— oo m=n— 2 

where we have recast recurrence relations in matrix form in Eq.(30), and defined matrix 
elements 

F - 



2D 



_ 2AB 
Ae + 

E - 

Condition that solutions of simultaneous equations in Eq.(30) exist is that deter- 
minant det(£') must be zero, that is 

det(E) = 0, (31) 
det|E-/ie/| = 0. (32) 

Secular equation of Eq.(32) is a characteristic equation that determines the exis- 
tence of periodic solutions of Eq.(29). Solutions R{^) could be functions with period 
47r/^. There exist four series of periodic functions according to period being odd or 
even. Eq.(27) has two confluently irregular singular points x = 0, oo. Behaviors of its 
solutions at event horizon Vh — 0(x = 0) depend upon that of R{^) at ^ — > ±ioo ( 
according to v being negative number or positive number ). 
V. CONNECTION BETWEEN THE RADIAL EQUATION IN NON- 
EXTREME CASE AND THAT IN EXTREME CASE 
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In this section, we illustrate that the radial equation of Eq.(27) in the extreme case 
is a confluent form of Eq.(13) in the non-extreme case, and give expression to the first 
thermodynamic law in the extreme KNBH case. 

After making substitutions of e = Me, ez = x, ez = Mx, B = eB, (0 < e < 1) in 
Eq.(13), we have 

r = M + ez = M(l + x),A = e\z'^ - 1) = M^x^ - e^), 

A = 2Mlj - qe, eB = MB = uj{2M^ - e^) - qeM - ma. 
then, Eq.(13) is equivalent to the following one in the non-extreme case: 

^[{x' - e')^^] + [M'ico' - fi')ix' - e') + 2M{Au; - M^,')x 

^ {Ax + Bf ^ _ ^2W2^2 _ g2^ _ ^q^j^^ _ x]R(x) = 0. (33) 

x^ — 

Eq.(33) has two regular singular points x — ±e, {z — ±1) which are located at 
exterior horizon and interior horizon ( Cauchy surface) r± — M±£ = M(l±e) 
respectively, along with another irregular singular point x = oo. After taking limits 
e — > 0, — — > x^, Eq.(33) in the non-extreme case tends to Eq.(27) in the extreme 
case. The latter has two confluently irregular singular points x — 0,oo. The irregular 
singular point x = which is located at event horizon r/j = M in the extreme case is 
just one to which two irregular singular points x — e and x — —e in the non-extreme 
case concur when e or £ — > 0. 

In the extreme KNBH case (M^ = a"^ + e^), surface gravity Kh = 0, event horizon 
r/j = M, reduced event horizon area Ah — M'^-\-a? — 2M^— e^, the first thermodynamic 
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law of extreme Kerr- Newman black hole is expressed as follows 

dM = VLhdJ + ^hde (34) 

where $/j = {erh)/Ah, and flh = a/Ah are electric potential, angular velocity at event 
horizon (r^ = M) respectively. 

VI. CONCLUSION 

In this paper, a sourceless charged massive scalar Klein-Gordon field equation has 
been separated into the angular and radial parts. The separated equations are all gen- 
eralized spin-weighted spheroidal wave equations. In the non-extreme case, we present 
general solutions in power series expansion and that of integral forms, as well as sev- 
eral special solutions with physical interest for the radial equation. These solutions 
can be orthonormalized to a set of complete functions. In addition, they have asymp- 
totic behaviors of plane waves at infinity. On the base of these orthogonal functions 
or polynomials, we can expand wave function of a complex scalar field to a quantized 
Klein-Gordon field on the Kerr-Newman background. In the extreme case, the radial 
equation can be reduced to modified Whittaker-Hill equation. In both cases, we obtain 
five-term recurrence relations between coefficients in power series expansions. 

At base of this work, the quantum conservation laws about Hawking process and 
probable generalization to black hole thermodynamic laws can be discussed further. It 
is anticipated that the separated parts of Dirac equation in the Kerr-Newman geometry 
could be reduced to the forms of generalized spheroidal wave equation. 
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APPENDIX 

In this appendix, we present three-term recurrence relations between coefficients 
in power series expansions around regular singular points {z — ±1) for spin- weighted 
spheroidal wave equation of Eq.(18), namely 

(1 - z^)W^{z) - 2[a +{b + l)z]W;,{z) + [c\z^ - 1) + 2dz + A„ - b]Wn = 0, 

{-l<z< 1) (Al) 

Eq.(Al) has two regular singular points z — 1 and z — —1, with indices p_ = 
0, —a — b and p+ — 0,a — b respectively. When c,d^O, Wn{z) must tend to Jacobi 
polynomials. 

Introducing a symbol e = ^1, we denote these two regular singular points z — ±1 — 
—e. Then, we make power series expansions around regular singular points z — —e 

respectively, where we have written them in an united manner: 

oo 
n=0 

Substituting the above regular solutions of Eq.(A2) into Eq.(Al), we obtain three- 
term recurrence relations between coefficients as follows 

(1 + 6 - ea)/i + [Ao + 2ac-b- 2e{bc + c + d)]/o = 0, 
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(n + l)(n + 1 + 6 - ea)fn+i + [A„ + 2ac-b- 2e{hc + c + d) 
-n(n +l + 2b + 4ec)]/„ + 4e(nc + bc + d)fn-i = 0. 
After defining coefficients, 

An = A„ + 2ac- 6 - 2e(6c + c + (i) -n(n + 1 + 26 + 4ec), 
Bn = (n+ l)(n + 1 + 6- ea), 
Cn = 4e[(n + 6)c + d], 

recurrence relations for first term and n-th term can be written as 
5o/i + ^o/o = 0, 

Bnfn+l + Anfn + Cnfn-1 = 0. (^3) 

Three-term recurrence relations of Eq.(A3) can handled by continued fraction 
method, ^'^ or by matrix method ( see J. W. Liu's paper in Ref.7 ) as we can ar- 
ray An,B„, Cn to make up a generalized Jacobi tridiagonal band matrix. Similar 
three-term recurrence relations can also be obtained by expansions in the light of Ja- 
cobi polynomials, but the coefficients An, Bn, Cn will be more complicated than those 
presented here. 

The second regular solutions around the same points can be easily obtained by 
Probenius's method, and we have not presented them here. Irregular solutions are 
connected with these regular ones by integrals similar to those in Eqs.(19) and (20). 

In order to make Wn{z) finite at z = ±1, Wn{z) must be truncated to be polynomi- 
als, then Wn{z) is orthonormahzed with eigenvalue and weight (1 — z)^'^^ {1 + zy~°' . 
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Hence we have 

Battle-Lemarie wavelet or Daubechies' compact support wavelets^^ can be used in 
numerical computation for matrix equation of Eq.(A3) and to prove convergence of 
polynomials Wn{z), but we don't pursue this goal here. 
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